Rotational Brownian motion on the sphere surface and rotational relaxation 
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The spatial components of the autocorrelation function of noninteracting dipoles are analytically 
obtained in terms of rotational Brownian motion on the surface of a unit sphere using multi-level 
jumping formalism based on Debye's rotational relaxation model, and the rotational relaxation 
functions are evaluated. 
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Rotational Brownian motion and rotational relaxation 
phenomena has received considerable interest in the lit- 
erature for a long time. Historically, the problem of 
isotropic rotational Brownian motion arose in connection 
with the work of Purcell on spin relaxation in systems 
with more than two spins in the same molecule. The dy- 
namics of rotational Brownian motion of a sphere around 
a given axis was discussed briefly by Einstein [l| in one of 
his early papers. Later, Debye |2| proceeded by extend- 
ing Einsteins treatment of the translational Brownian 
motion to rotational Brownian motion of noninteracting 
permanent dipoles in the presence of an external time- 
varying field and solving the appropriate Fokker-Planck 
equation. Here that equation is the Smoluchowski equa- 
tion, which is an approximate Fokker-Planck equation 
in the space of angular coordinates for the distribution 
function of the orientations of the dipoles on the surface 
of the unit sphere when the influence of the inertia of 
the molecules on the relaxation process is ignored. Fol- 
lowing years many generalizations of Debye model have 
been extensively studied 0, S IE IE S H, S HO, How- 
ever, Debye theory cannot explain relaxation phenomena 
in strongly interacting systems. The relaxation behav- 
ior in strongly interacting systems such as amorphous 
polymers, spin glasses, glass forming liquids, etc., devi- 
ates considerably from the exponential (Debye) pattern 
$ ~ exp [— t/r] and is characterized by a broad distri- 
bution of relaxation times. The relaxation in strongly 
interacting systems obey to the stretched exponential re- 
laxation function experimentally $ ~ exp [— (t/r) a ] with 
< a < 1, which often referred to as the Kohlrausch- 
William- Watts (KWW) function [IE IH 

In this Letter our aim is to obtain rotational relax- 
ation function in terms of autocorrelation function of ro- 
tational Brownian motion using multi-level jumping pro- 
cess (MJP) formalism based on historical Debye rota- 
tional relaxation scenario. The MJP is the generalization 
of the two-level jumping process (TJP) formalism. The 
TJP and MJP have been used to obtain the relaxation 
function 

HI 13 13 113 Recently, the non-exponential 
relaxation function in complex systems has been obtained 
using MJP formalism [Isj . 



Autocorrelation function <p is simply given by 
0(*) = </z(O)/z(*)) 



(1) 



where /i(0) and fi(t) denote values of fluctuation variable 
fi, i.e. dipole, at time zero and the same variable at a later 
time t, respectively. Physically, it measures the time over 
which the variable fi retains its own memory until this 
memory is averaged out by statistical randomness. 

If it is assumed that the external field has been applied 
in the z-direction, the spatial components of the autocor- 
relation function of dipoles can be written in terms of the 
Brownian motion on the surface of the unit sphere in a 
simple form 

(j) x (t) = (fi x (0) fi x (t)) = (fi sin 9 cos <p) (2a) 



<\> z (t) = (fi y (0) fi y (<)) = (fio sin 9 sin ip) (2b) 



<f> z (t) = (fi z (0) fig (t)) = (no cos ( 



(2c) 
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where fiQ denotes magnitude of the dipole (/io = 1), 
and ip are stochastic variables which define the spatial 
orientation of a dipole, in other words, position of the 
Brownian particle on the sphere surface. For simplicity, 
it can be considered that both of 9 and if are discrete 
stochastic variables. The 9 takes random values in be- 
tween and 7T, and on the other hand, <p takes random 
values in between and 2ir, where 7r's value is defined 
in terms of degree. Hence, stochastic variable fi can take 
randomly integer or non-integer values depending on 9 
and ip. In a such process, fluctuation variable fi jumps 
from one value of fi to another stochastically. As it can 
be seen in Eqs. (2) only z-component of the autocorrela- 
tion function depends on 9, however, x and y components 
depend on two variables i.e. 9 and if. We will see below 
that all components of the autocorrelation function are 
carried out using MJP after a little bit of algebra. 

All components of the autocorrelation function in 
Eqs. (2) can be derived by using stochastic formalism. 
For example, the z-component 4> z (t) of the autocorrela- 
tion function is given by 

^ z (*) = X! p VzoP (VzO | fig, t) fl z (3) 
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where the quantities [i z § and fi z are the values of stochas- 
tic variable \i z at times and t, respectively. Having 
started at the initial state /j, z q with the statistical weight 
p(Mzo); and P {fJ-zO I fJ-z,t) measures the probability of 
propagation from /jl z q to fj, z in time t, which is known as 
conditional probability. 

The conditional probability P (fJ, z o \ fJ> z ,t) can be given 
by the matrix representation associate a stochastic state 
\fi z ). In the operator formalism (See Refs. 0,E|) the 
conditional probability P (fJ- z o I > t) is simply written 
with the matrix elements of operators as 



P {VzO \ Vz,t) = (Hz\P (t) \Hzo) 



(4) 



in the matrix form where P(t) is the conditional probabil- 
ity operator, and |.„) indicates vector space of stochastic 
states. Note that the conditional probability satisfy 

P{Li z o | p„0) = (fi z \P{0) \fizo) =*(/**- Mzo) (5) 
for t = 0. 

Similarly, the fluctuation variable /i z can be also given 
by matrix form, the eigenvalues of which correspond to 
stochastic variables /io, H>i,- respectively. The ma- 
trix representation the non-interacting fluctuation vari- 
able n z can be represented associated with operator Z as 



Z 



/Mo 
/l! 
fl 2 



V 



\ 






(G) 



Eigenvalues of operator Z are given by \i z —cos9 where 

•- 0,1,2 -• 

Thus, the autocorrelation function in Eq. (3) can be 
transformed to the matrix representation 

<t ) z(t)= Pfazo) Z \^zo) (fJ-z\P(t) \Z Z0 ) (fJ, z \ Z\H Z 

(7) 

using above definitions. 

As yet, the autocorrelation function (3) has been given 
by matrix representation as seen in Eq. (7). The Eq. (7) 
can be solved easily using matrix algebra and with more 
convenient basis vector notation if the operator form of 
the conditional probability is known. 

The operator form of the conditional probability can be 
evaluated from master equation. The conditional proba- 
bility P (fJ-za I l^zjt) can be given associated with master 
equation 



d_ 



P (flzO \tiz,t) = Y, W I P (M*o I M*,*0 (8) 
Mi 



where W (fi' z | fj, z ) is the jump rate from fj, z to [i' z of the 
Brownian particle i.e. dipole. 



The jumping rate W (fi' z \ fi z ) can be given by matrix 
form 



W(fx' z \»*) = (^\W\» z ) 



(9) 



where W is the jumping rate operator. Thus Eq. (8) can 
be rewritten completely by a matrix form as 

§- t (Vz\P (t) |MzO> = X>* I ^ W, I P W lMzO> 

(10) 

using Eq. (4) and (9). 

The stochastic states taken to form an orthonormal set 
which provide closure property is as 



(ii) 



for discrete variables. Using closure property, Eq. (10) is 
reduced to 



dP(t) ^ 

(» z \^T 1 \Li z0 ) = {» z \WP(t)\ f i z0 ) 

This equation can be simplified as 
d 



dt 



p (t) = wp (t) 



As a result, the solution of Eq. (13) is 
P (t) = exp (wtj 



(12) 



(13) 



(14) 



where W is the jump matrix and can be written in terms 
of collision and unit matrix elements as 



W = A J- 1 



(15) 



where A is defined as the relaxation rate, J is the collision 
matrix and 1 is a unit matrix. The relaxation rate A for 
7r-level jumping is given by 

A = irw 



(16) 



and collision matrix J is presented by 



/ 1/tt 1/tt 
1/tt 1/tt 



J = 



1/tt 



1/tt/ 



(17) 



V 1/tt 1/tt 

where w is the jump rate from one value of [i z to another. 

It is helpful to write down the explicit form of the jump 
matrix W in terms of w. If Eqs. (16) and (17) are inserted 
into Eq. (15), the jump matrix W can be constructed as 



W 



/ — (tt — l)w 
w 



w 

-(•7T — l)w 



w 



w 
w 



\ 



-(w — l)w / 



(18) 
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The usefulness of the decomposition in Eq. (15) comes 
out from the fact that the J matrix has a very simple 
property 



J — J • J — J • ' ' ' j J ' — J 



(19) 



for any integer k > 0. This property allows us to imme- 
diately construct the conditional probability matrix as 



P (t) = cxp j^A (J - 1J tj (20) 
or, alternatively it can be given in the matrix form 

n 



( r] + e 



-At 



P(t) 



r) + e 



-At 



V v 

where 77 is defined by 



. . . 77 


\ 


. . . 77 




. . . rj + e' 


-At j 



(21) 



r?=-(l 

7T 



-At 



) • 



(22) 



The conditional probability operator P (t) can be ex- 
pressed in useful form as 



P (t) = cxp (-At) 1 - J + J exp (At) 



(23) 



using direct power series expansion and relation which is 
given by Eq. (19). 

On the other hand, it is convenient to introduce the 
more general notation for basis vectors to solve Eq. (7). 
Therefore, let us associate a stochastic state \8) (9 = 
0, 1, 2, 7r for multi-level jumping process) instead |/_i z ) 
of 7r values 



|0> 








,11} = 



1 





w 



,...,k) 








Vi/ 



(24) 



The closure property for new stochastic states \0) (and 
similarly \9')) is presented as 



Ei^) <*i = 1 



(25) 



in a discrete form. 

Now, Eq. (7) can be rewritten using new basis vectors 

as 

<p z (t) =Y,Pe (0\ Z \9) (6'\ P (t) \9) (9'\ Z \9') (26) 



with apriori occupation probability 



1 



P0 = -, fl = 0,l,2,...,7T . 
7T 



(27) 



The expectation value of conditional probability ma- 
trix P (t) in Eq. (26) is obtained as 



(e\P(t)\e') = - + 

7T 



- ) exp (-At) (28) 

7T , 



using Eq. (23) and new basis vectors in Eq. (24), and on 
the other hand, the expectation value of Z is given by 



{6\Z\e') = Z e 5ee> 



(29) 



with 9, 9' = 0, 1, 2, 7r, where Zg are the allowed values 
of the stochastic variables. 

If Eqs. (27)-(29) are inserted into Eq. (26), the autocor- 
relation function (26) can be reduced to a simple form 

0, (t) = {^f + {{A) - {^f) cxp (-At) (30) 

where the deterministic quantities arc weighted averages 
over the available states of the corresponding variable. 
The average values of (/j z ) and (/j z ) in the stationary 
state are given by 



<Mz) = I> {e\Z\6) = lY j Z e 



(31) 



0=0 



e=o 



and 



(£) = j2p n m 2 \9) = lY.zj (32) 



9=0 



8=0 



respectively. The diagonal elements of Z matrix (6) are 
given as fi z — cos#, i.e., Zg = cos9. Thus, it is expected 
that (fj, z ) = and = 1/2. 

Finally, z-component of the autocorrelation function 
for the Brownian motion on the spherical surface can be 
obtained simply as 



(t) = - cxp [-At] 



(33) 



As it can be seen in Eq. (33) z-component of the ro- 
tational autocorrelation function has exponential form 
which indicates that the same component of rotational 
relaxation must has exponential form. 

The relaxation function for z-component can be imme- 
diately obtained using linear response theory 



* (t) = * (t = 00) - * (t) 



(34) 



where ^ (t) is well known response function which can be 
associated to the famous response-relaxation function as 



*(*) 



1 

It 



((^)-(mz (o)Mt)» 



(35) 
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where k is the Boltzmann constant, and T is the tem- 
perature. If Eq. (35) is inserted into Eq. (34), as a final 
result, the rotational relaxation function can be obtained 
as 



2kT 



exp [—At] 



(36) 



We have shown that the z-component of the autocor- 
relation function and relaxation function can be obtained 
using the multi-level jumping formalism. Similarly, the 
other components of the autocorrelation function, that 
is, cj) x (t) and 4> y (t) can be also derived using this tech- 
nique. However, to obtain the components <f) x and <p y , 
we should eliminate the ip in <j) x and <f> y . 

At this point, firstly, it can be suggested that com- 
ponents of the autocorrelation function in Eq. (2) must 
satisfy the relation 

<0l(*)> + (#(*)> + = ! (37) 

which is caused from a simple trigonometric identity as 
(sin 2 6 cos 2 <p) + (sin 2 9 sin 2 <p) + (cos 2 6) = 1 . (38) 
Therefore, the components <p x and <p y can be written as 
(^)+^(*)> = l-(cos 2 0) . (39) 

On the other hand, it is expected that the component 
4>y it) should be equal to (f> x (t). Thus, <f> x (t) is written 
as 



X {t) = (<f> x (t)) = -j={awl 



(40) 



in simple form. 

As it can be seen in Eq. (40), using mathematical trick, 
we have got rid of ip in the components <j) x (t) and 4> y (t) . 



Now, to calculate the component cj> x (t) we need only val- 
ues (n x ) and (/i 2 ). The averages can be given as 



7T 7T 
L < n £ / 



(41) 



9=0 



and 



7T 7T 

{nl)=Y.P^ e \X 2 \ e ) = -Y. X e ( 42 ) 



9=0 



8=0 



respectively. The matrix X has the same form with (6). 
But elements of X matrix are equal to \i x = sin0, i.e., 
Xg = sin 8. Therefore, required averages can be carried 
out as (fi x ) = \/2/tt and (// 2 ) = 1/4. Hence, the re- 
laxation function for x and y components are given by 



$, (t) = % (t) = -L Q - cxp [_ At ] . (43) 

The Eq. (43) clearly show that the components (j> x (t) and 
cj) y (t) of rotational relaxation function have also exponen- 
tial form. 

In this Letter, the spatial components of the autocorre- 
lation function of noninteracting dipoles have been ana- 
lytically obtained in terms of rotational Brownian motion 
on the surface of a unit sphere using multi-level jumping 
formalism based on Debye's rotational relaxation model, 
and the rotational relaxation functions have been evalu- 
ated. As it is expected, all components of the rotational 
relaxation function have exponential form when the in- 
fluence of the inertia of the dipoles is ignored. 
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